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We forecast the ability of cosmic microwave background (CMB) temperature and polarization 
datasets to constrain theories of eternal inflation using cosmic bubble collisions. Using the Fisher 
matrix formalism, we determine both the overall detectability of bubble collisions and the constraints 
achievable on the fundamental parameters describing the underlying theory. The CMB signatures 
considered are based on state-of-the-art numerical relativistic simulations of the bubble collision 
spacetime, evolved using the full temperature and polarization transfer functions. Comparing a the¬ 
oretical cosmic-variance-limited experiment to the WMAP and Planck satellites, we find that there 
is no improvement to be gained from future temperature data, that adding polarization improves 
detectability by approximately 30%, and that cosmic-variance-limited polarization data offer only 
marginal improvements over Planck. The fundamental parameter constraints achievable depend on 
the precise values of the tensor-to-scalar ratio and energy density in (negative) spatial curvature. 

For a tensor-to-scalar ratio of 0.1 and spatial curvature at the level of 10~^, using cosmic-variance- 
limited data it is possible to measure the width of the potential barrier separating the inflating false 
vacuum from the true vacuum down to Mpi/500, and the initial proper distance between colliding 
bubbles to a factor 7r/2 of the false vacuum horizon size (at three sigma). We conclude that very 
near-future data will have the final word on bubble collisions in the CMB. 


I. INTRODUCTION 

Imminent results from the Planck satellite [T] will contain nearly all of the large-scale cosmological information 
encoded in the temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. 
This makes the present a highly opportune time to determine which analyses can maximize the scientific return from 
this dataset. One exciting opportunity is the potential observation of relics left from events that occurred in the very 
early Universe. In this paper, we forecast the ability of a cosmic-variance-limited CMB dataset to constrain one such 
class of early Universe events: cosmic bubble collisions in eternal inflation. 

Inflation, a hypothesized epoch of accelerated expansion, has become an essential component of the standard 
cosmological model mg. As a side-efFect, inflation can in many cases give rise to an eternally inflating multiverse 
[i[7]. In this scenario, a high-energy inflating phase is exited locally, inside bubbles, but not globally. In eternal 
inflation, the rate of bubble formation is outpaced by the accelerated expansion of the high-energy inflating phase, 
preventing the percolation of bubbles and leaving an increasingly large volume in which inflation and bubble formation 
continues. An important test of this scenario is the observation of the wreckage left from collisions between our own 
bubble and others [8]. Determining the outcome of collisions and the probability of observing them has been the 
subject of a substantial body of work [9l[42]. Observational constraints on this scenario have already been placed 
using temperature data from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite [431148] . 

There have been a number of important developments since the early work on constraining bubble collisions with 
observations. Importantly, a direct link has been made between the scalar field Lagrangian and cosmological observ¬ 
ables [24l |28] . For models containing a single scalar field with a canonical kinetic term, the properties of the eternally 
inflating Universe are determined entirely by the potential of the scalar field. An example is shown in Fig. this 
model allows two types of bubbles. Because a bubble collision gives rise to an inhomogeneous Universe, observers at 
different locations will have access to different parts of the collision spacetime. For observers near the causal boundary 
of the collision there is a simple analytic template for the comoving curvature perturbation caused by a single bubble 
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collision in single-scalar-field models pS], given by 






r-obs 


Mpi 


(1 cos ^Xgep) 


X — Xc 


^Is 


+ 


-^coll ^coll 


-obs 


Rohs 


cos 


2 (a; - XcY 


( 1 ) 


where is the curvature inside the bubble containing the observers, and are the inflationary tensor-to- 
scalar ratios inside the observation and collision bubbles, and are the Hubble scales during inflation inside 
the observation and collision bubbles, width of the potential barrier separating the inflating false vacuum 

from the true vacuum in the collision bubble, Mpi is the Planck mass, 0 < Axgep < tt is the distance between the 
colliding bubbles (measured in terms of the false-vacuum Hubble parameter) in the centre of mass frame, Xc is the 
comoving position of the causal boundary of the collision, and x\s is the comoving distance to the surface of last 
scattering. The parameter Axgep varies from collision to collision, and Xc depends on the position of the observer; 
all other parameters are fixed by the scalar field Lagrangian as shown in Fig. Eq. [^is valid in the limit where the 
slow-roll approximation holds in the future of the collision inside both the observation and collision bubbles. Outside 
this limit, numerical simulations are necessary to accurately determine the template. All cases considered below are 
well described by Eq. 



FIG. 1: An example of a potential giving rise to eternal inflation with two possible types of bubbles. The components of the 
potential determining the various parameters in Eq. [^are labeled. 


The collision spacetime possesses SO(2,1) symmetry, which translates into an approximate planar symmetry in the 
neighbourhood of an observer in the limit <C 1. This geometry is shown in Eig. The planar symmetry of the 
collision and the existence of a causal boundary separating the bubble into regions that are or are not affected by the 
collision imply that the effects of the collision are confined to a disc of angular radius Oc on the sky of an observer. 
The mapping between the position of the causal boundary Xc (in Mpc) and the angular scale of the collision for a 
frame in which the observer is at the origin of coordinates is shown in the right panel of Eig. 

In this paper, we forecast the ability of a suite of CMB experiments, producing both temperature and polarization 
data, to detect bubble collisions of the form given in Eq.[^ We forecast constraints on the overall detectability of bubble 
collisions as well as on the fundamental parameters underlying the collision model. Although constraints [43H48] and 
forecasts [22] exist for bubble collisions in the CMB, so far they have assumed a phenomenological and incomplete 
template. This paper is an important step forward as we forecast constraints directly on the scalar-field potential 
that underlies eternal inflation. As the relevant datasets become available, searches along the lines of Refs. |43H48| 
can be implemented using the complete collision template. 

The paper is structured as follows. In Sec. |n| we describe how the curvature perturbation template is evolved to 
the CMB observables. Sec. [TTT| describes our forecasting technique, which is based on the Eisher information matrix. 
In Secs. [IV| and [V| we present our forecasts and compare them to existing constraints; we conclude in Sec. [VT 


^ In generating this plot, and throughout the paper, we have assumed the 2013 P/anc/c-hWP-hhighL-hBAO best fit cosmology m- 
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FIG. 2: In the left panel, we show the geometry of a bubble collision on the surface of last scattering. Shaded regions are 
affected by the collision, while unshaded regions are not. The observer is located at the origin, and the circle represents their 
past light cone. In the right panel, we show the mapping between the observed angular radius of the collision Oc in degrees and 
the comoving position of the collision boundary Xc in Mpc. 


II. EVOLVING THE PRIMORDIAL COMOVING CURVATURE PERTURBATION 

Given the analytical parameterization of the curvature perturbations in Eq. we now describe the formalism 
adopted to predict bubble signatures imprinted onto the CMB. The spherical harmonic coefficients 6^^ of the template 
are related to the curvature perturbations 7lirn{k) via (e.g., Ref. [50] ) 

bfm = ^ , ( 2 ) 

where Af-{k) is the radiation transfer function in temperature and polarization {X = {T, E}) in momentum space.^ 
For an efficient solution of the problem, it is more convenient to phrase the equation in terms of comoving distances, 

= J dr af (r) Tlem (r), (3) 

where we have introduced the real-space transform of the radiation transfer function, 

J ^eik) je{kr), (4) 

where the ji{kr) are spherical Bessel functions of order £. 

To evaluate Eq. numerically, we first obtain the radiation transfer function using a modified version of the 
Boltzmann integrator CAMB [51]. In doing so, we assume the 2013 P/ancA:+WP+highL+BAO best-fit cosmological 
parameters m, adopting the limiting case of a fiat cosmology. Given our parametrization of as a function of 

we then compute the radial integral on 90 nodes for all multipole moments \m\ < 2500. We use the approach 
introduced in Ref. [52| to optimize node positions and quadrature weights to minimize the error associated with the 
numerical integration. Resolution studies have proven results to be stable. 

In Fig.[^ we visualize the outlined procedure for the linear term IZ oc x—x^^ with the distance to the causal boundary 
of the collision x^ set to 10 Gpc. We plot the curvature potential on a three-dimensional spherical grid enclosing 
the observable Universe in comoving coordinates (left), and juxtapose the derived templates for the temperature and 
polarization intensity (defined by P = + U‘^) signals (right). Whereas the former is a relatively smooth function 

with a well defined boundary, the latter shows a more extended signature and a comparatively sharp feature at the 


^ In our sign convention, the CMB temperature anisotropies in the Sachs-Wolfe approximation are related to the comoving curvature 
perturbation by 5T/T = IZ/b. 
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FIG. 3: The induced signal in the three-dimensional curvature perturbation {left panels shown on shells centered on the observer 
out to the last scattering surface) results in radially symmetric features in the CMB. Right panel: with the direction of the 
bubble collision in the center of the plot, we show Mollweide projections of the corresponding templates in temperature {upper 
half) and polarization {lower half) in dimensionless units. 


location of the observed angular radius of the collision. This feature is a result of the discontinuous first derivative of 
the curvature perturbation at x = Xc [22], and is consequently absent in the quadratic term, where IZ (x {x — Xc)^. 

In practice, we need only transform the two position-dependent terms {x — x^) and {x — Xc)‘^ individually to build 
any template desired (at a given scale). We normalize the curvature perturbations so that they are fixed to 0 at x = Xc 
and 1 at X = xis, and denote their evolved forms Lf^{xc,uj) and Qf^{xc^oo)^ respectively, where uj = {O^cf) are the 
angular coordinates of the centre of the signature on the sphere. To minimize the number of parameter degeneracies 
we face, we restrict our analysis to collisions between identical bubbles, setting and In this 

case, the spherical harmonic coefficients of a bubble collision in the CMB temperature or polarization can be written 
as 
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( 5 ) 
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where X = {T, E} and the constants 


A = 



B = 20f ® 


( 7 ) 


are assumed to have been constrained by other data. When presenting constraints on fundamental parameters, we 
assume the hypothetical situation where primordial gravitational waves and negative spatial curvature have been 
detected at a level of = 0.1 and = 10“^, respectively, yielding the numerical values A 0.179 and 

B = 2 X 10 Throughout our analysis, and without loss of generality, we choose to centre the bubble collisions on 
the North Pole. This, combined with the azimuthal symmetry of the bubble collision signature, means that all 
with m 7 ^ 0 are zero. 
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Experiment 

V / GHz 

aF / \jK 

oF! \jK 

Qv / arcmin 

Planck 

143 

4.1 

7.8 

7.1 

217 

8.9 

18.2 

5.0 

WMAP 

94 

35.3 

- 

13.2 


TABLE I: The experimental characteristics assumed in this work. 


III. FISHER MATRIX ANALYSIS 


Having obtained the signature of bubble collisions in the CMB, we now describe the method whereby we forecast 
constraints on the underlying parameters. We assume that the data — in this case, the observed spherical harmonic 
coefficients where X = {T^E} — consist of the stochastic Gaussian CMB the beam-deconvolved in¬ 

strumental white noise (n^), and a deterministic bubble collision (6^). Under these assumptions the moments of 
the (Gaussian) data are simple to define: the noise and CMB anisotropies do not contribute to the mean, which is 
determined entirely by the collision to be 


(d) = n 



(8) 


and the deterministic bubble collision signature does not contribute to the covariance, which is defined purely by the 
CMB and noise power spectra to be 


(dd*) - (d)(d)‘ = c = 


CJT + ^rr 
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( 9 ) 


We generate the CMB power spectra Ci using CAMB, fixing the cosmological parameters to their 2013 
P/ancA;+WP+highL+BAO best-fit values [49] 

We produce forecasts for three experimental configurations: a theoretical cosmic-variance-limited experiment with 
= 0, a forecast of the Planck satellite’s full mission, and WMAP, included to compare our forecasts against existing 
results. Following Ref. [53|, we define the Planck forecast to correspond to 30 months of observations, retaining only 
two CMB channels (143 and 217 GHz) to account for foreground subtraction. The beam-deconvolved noise power 
spectrum for this mission is 






81n2 



( 10 ) 


where is the full-width at half-maximum of each channel’s (Gaussian) beam in radians and is the noise 

variance per beam-sized patch in temperature or polarization. Our WMAP experiment corresponds to the seven- 
year, foreground-reduced W-band temperature map [54]. The precise values of the experimental parameters are 
presented in Table [TJ Note that we do not include the effects of partial sky coverage in our analysis, and the forecasts 
we present therefore apply to signatures that are not significantly masked.^ 

The scope of this work is restricted to investigating the parameters describing the bubble collisions (6^^), so the 
covariance is independent of the parameters of interest. In this setting, the Fisher matrix is given by [55] 
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Assuming the likelihood C is Gaussian in the bubble parameters, we can associate the inverse of the Fisher matrix 
with the covariance of the parameters 




( 12 ) 


^ The mode-count reduction due to masking can be approximated by inserting a factor of in the covariance matrix (Eq. 9l. The 

_ 1/2 

net effect is to boost all parameter uncertainties by a factor of . Note that the relevant sky fraction depends on the size of the 

collision and its position with respect to the galactic plane and other masked sources. 














6 


The diagonal of the inverse Fisher matrix corresponds to the variance on each parameter after having marginalized 
over all others; the inverse of each diagonal element of the Fisher matrix is the variance on each parameter conditional 
on all other parameters taking their maximum-likelihood values. 


A. Parameterizations and Detectability 


Our goals in this work are to forecast the detectability of bubble collision signatures as well as our ability to 
constrain the fundamental parameters describing them. Both of these goals can be achieved using the Fisher matrix 
formalism described above, simply by re-parameterizing the collision signature. The simplest starting point is to cast 

C / coll 

the analysis in terms of the Lagrangian parameters {0 = Axgep, The resulting Fisher matrix encodes 

our ability to measure the fundamental parameters describing the bubble collision space-time, but does not yield a 
simple concept of detectability. To define our detectability criteria, we re-cast the templates in terms of observable 
amplitudes, and require that these amplitudes can be distinguished from zero with some threshold significance. 

With this concept in mind, we employ two additional parameterizations in this work, based on observable rather 
than fundamental quantities. Specifically, we re-cast the template in terms of the individual amplitudes of the linear 
and quadratic terms {O' = {71q , Xc^oo}), or the total amplitude and the fraction of the template contributed by 

the linear term {0" = {TIq, The various amplitudes are defined by 


bfm = ^0W) + n^Qf^{Xc,U) 
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and their parameters are related via 
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By determining at what point the overall (T^o): linear (T^q )? quadratic {IZq) amplitudes become distinguishable 
from zero at a given significance, we are able to state when mixed, pure-linear, and pure-quadratic bubble collision 
signatures can be detected. 

To obtain constraints on the full range of parameters we need only calculate the Fisher matrix in terms of one set 
and transform the resulting inverse matrix using 




Y. 
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kl 
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(18) 


Fundamental to the Fisher matrix formalism is the assumption that the likelihood is Gaussian in the parameters of 
interest. This guarantees that the initial parameterization choice has no bearing on the Fisher matrices produced for 
each other parameterization; it also means that the resulting forecasts are more accurate for parameterizations which 
better satisfy this assumption. As the likelihood is bivariate Gaussian in the linear and quadratic amplitudes, we 
expect these to be most accurately predicted, whereas the error for the total amplitude, for example, will be larger. 


B. Implementation 

As our aim is to forecast our ability to detect and constrain the bubble collision model as a function of the Lagrangian 
parameters, we frame the Fisher matrix analysis in terms of this parameterization, with one small change: in order 
to linearize the dependence of the template on the parameters (and hence make the likelihood as close as possible to 
Gaussian in the parameters), we use 1 — cos Axgep instead of Axgep- From Eq. we see that the derivatives with 

^ ,coll ^ 

respect to and 1 — cos Axgep can be computed analytically. We calculate the derivatives with respect to x^ 

numerically, using two-sided finite differences with a step size of ±1 Mpc. We have checked that the derivatives are 
stable over a range of adjacent step sizes. The azimuthal symmetry of the collision signatures and our chosen central 

e , coll 

position guarantee that the derivatives with respect to , 1 — cos Axgep and x^ are non-zero for m = 0 only. 
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Though the templates in general depend on the angular position of the centre of the bubble collision, for our 
particular choice of centre {uj = (0,0)) we can neglect both coordinates from the Fisher matrix analysis. Firstly, the 
derivative with respect to the longitude 0 is zero due to the azimuthal symmetry of the signatures. Secondly, the 
derivative with respect to the colatitude 0 — formed by rotating the templates to 0) and calculating two-sided 

finite differences — is non-zero for odd m only. As a result, all off-diagonal entries in the Fisher matrix involving 
0 are zero, and 0 is therefore entirely uncorrelated with the parameters of interest. The Fisher matrix we calculate 
therefore reduces to 


dbj^ ^eE 

E dOi do n ^ 

-^- 


dbfo db] 
dOn dOn 


10 


dbfo db^ 
dOi dOn 


10 




TE 

£ 


(jTT(jEE _ ^qTE^^2 


9bfo dbf^ ^tT 

dOi dOj ^£ 


(19) 


where we have absorbed the noise power spectra into Cf^. 

We evaluate the Fisher matrix on a 50x50x13 grid of parameter values, with linearly spaced samples in the ranges 

e i coll 

0.0005 < < 0.014, 0.01 < 1 — cos Axsep < 2.0 and 1000 < Xc < 13000 Mpc. At each sampled point in param¬ 

eter space we report the fractional marginalized uncertainty on each of the parameters of interest: the Lagrangian 

parameters 1 — cos Axgep and Xc, and the observable amplitudes 7Zq and TZq. 


C. Checks 


Although it is prohibitively slow to perform a brute-force likelihood analysis on this grid, we can evaluate the 
likelihood on a smaller grid to check our Fisher matrix outputs for the Lagrangian parameters. The likelihood we 
need to calculate is the probability of obtaining the data (which we assume contain a bubble, 6^^, with some fiducial 
parameter values) assuming our model (that the data contain a bubble, with a sampled set of parameters). 
To connect with the Fisher matrix analysis, we have to marginalize the logarithm of this likelihood over all CMB 
realizations. In the setting described above, this quantity takes the following form 
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where terms independent of the bubble collision parameters have been discarded. 

e / coll 

We evaluate this expression on a 50x50 grid of and 1 — cos Axgep, selecting Xc = 10 Gpc as a fiducial bubble 
size. The one-sigma conditional errors on each parameter are determined by finding the points at which the change in 

r i coll 

is one, linearly interpolating between grid points for extra accuracy. For we find sub-percent-level agreement 

between the errors derived via the Fisher matrix and the brute-force likelihood evaluation: this is to be expected, 

r i coll 

as the likelihood is Gaussian in -^y- For 1 — cos Axgep, in which the likelihood is not Gaussian, we find the Fisher 
matrix overestimates the uncertainty by factors of a few percent (when the error contours are well resolved by the 
grid) to factors of ^ 2 (where interpolation errors are large). In the regime where interpolating errors are small, we 
conclude that the Lagrangian parameter constraints presented below are accurate at the percent level despite the 
non-Gaussian nature of the likelihood for 1 — cos Axgep. The likelihood for IZq and 71q is Gaussian, and therefore we 
expect the Fisher matrix constraints to reproduce results obtained from the exact likelihood. 


IV. RESULTS 

Illustrative examples of the parameter uncertainties achievable for the cosmic-variance-limited experiment are plot¬ 
ted in Fig. for signatures with Xc = 1 Gpc (left) and 13 Gpc (right). Each panel contains coloured contours of 
constant fractional uncertainty in a Lagrangian parameter. Superimposed are shaded regions in which certain sets of 
the amplitudes are detectable at the three-sigma level. More precisely: 

1. in dark grey regions, none of the amplitudes are deemed detectable; 

2. in mid-grey regions, the total amplitude (T^o) is detectable, but the linear (T^q ) quadratic amplitudes 

are not; 












3. in light grey ^ hatched regions, the total and linear amplitudes are detectable; 

4. in light grey ^ hatched regions, the total and quadratic amplitudes are detectable; and 

5. in white regions, ah amplitudes are detectable. 

A brief comment on theoretical priors is appropriate at this stage. Assuming that bubble nucleation is a random 
stochastic process yields predictions for the theoretical priors over Xc and Axgep- The prior over Xc is uniform [l5l|20], 
and the prior over Ax^ep is proportional to sin^ Ax^ep na . Without a better microphysical theory of the potentials 

e / coll 

underlying eternal inflation, the prior over is unknown; however, a reasonable assumption would be to employ 

a logarithmic prior. In the signal-dominated regime with which we are primarily concerned, these priors do not 
signihcantly change the conclusions presented below. 


A. Template Detectability 


In general, the detectability boundaries (solid black contours in Fig. follow lines of constant amplitude 
or Tip). Note that the detectability boundary for the quadratic amplitude is missing from the Xc = 13 Gpc 
plot (and, indeed, for ah Xc ^ 5 Gpc): there is no combination of the Lagrangian parameters at these scales for 
which a pure-quadratic bubble collision signature is detectable. The reason for this is clear from the dehnition of the 
quadratic amplitude (Eq. 16). As Xc tends to xis, the term completely overwhelms the 1 — cos Axgep term. 


strongly suppressing TZq. Where it is present, the detectability boundary for the quadratic amplitude is, as expected. 


a function of 1 — cos Axsen only: for a given scale, IZq depends only on 1 


^sep 

/coll 


cos Axse^ This relation breaks down for 


C 7 coll I 

very small values of (a factor of ^ 10 smaller than the range plotted in Fig. 4), at which point a degeneracy 

opens up between 1Zq and Xc- Here, the detectability boundary turns upwards, becoming a strong function of . 
The sharper features of the linear template allow Xc to be better determined. 

Though IZq also contains a factor of (Eq. |l5k, we can always hnd a value of 

We 


it detectable. Unlike the quadratic amplitude, the *rmear amplitude is therefore detectabl^ in some portion of the 


large enough to make 




—(1 — cos Axsep) plane for ah values of Xc considered, albeit at higher and higher as Xc tends to xi^. As 

15) and Xc is well determined in the 
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the product of 


5^ 

Mpi 


coll 
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cosAxsep appears in the linear amplitude (Eq. 


parameter range considered here, the IZp detectability boundary is given by (1 — cosAxgep) oc const. 

To a hrst approximation, the total amplitude detectability boundary can be thought of as interpolating between 
the pure-linear and pure-quadratic cases. In fact, the linear and total amplitude boundaries cross over for very small 

shows that we should expect IZq and / to become completely degenerate. 


IZq . In this regime, examination of Eq. 
as we are trying to hx both of their values using only one number (Tip). The f parameterization is not a good 
model for regimes in which one of the templates is negligible. 
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B. Lagrangian Parameter Constraints 

We now turn our focus to the Lagrangian parameter constraints (the coloured contours in Fig.[^. Recall that these 
constraints depend on the detection of spatial curvature and primordial tensors, where we have assumed = 0.1 
and Of" = 10-^. 

First, we note that Xc (top) is the best constrained of the Lagrangian parameters, with percent-level errors achievable 
across much of the parameter space considered. Note that the shapes of the constant-Xc-error contours mirror that of 
the linear amplitude decision boundary: the sharp features in the linear template are critical in determining the size 
of the signature. 

r i coll 

The performance for the other two Lagrangian parameters, and 1 — cos Axgep, is much poorer, and highly 

dependent on the size of signature present. We only obtain direct constraints on 1 — cos Axgep through the quadratic 

r j coll 

template, with the linear template constraining the product [1 — cos Axgep]- This correlation structure is evident 
in the shapes of the individual contours. For signatures covering a large fraction of the sky we are able to measure 

e j coll 

both the linear and quadratic amplitudes, and hence constrain and 1 — cos Axgep with accuracies of up to ^ 10%. 
The smallest values of and Axgep measurable at three-sigma are roughly Mpi/500 and 7r/2, respectively. For 

smaller signatures, our inability to accurately measure the amplitude of the quadratic template, compounded by the 
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cr{xc) cr(xc) 




cr(S(/>o/Mpi) 








FIG. 4: Contour plots of the fractional marginalized uncertainty on the Lagrangian parameters {dashed, red to orange) for 
signatures with Xc = 1 {left eolumn) and 13 Gpc {right eolumn). Regions of parameter space in which different components 
are distinguishable from zero at 3cr are indicated by shading. In dark grey regions, none of the amplitudes are detectable; in 
mid grey regions, IZo is detectable but IZq and IZq are not; in light grey ^ hatched regions, IZq is detectable in addition to 
7^0; in light grey ^ hatched regions, TIq is detectable in addition to TIq. In white regions all amplitudes are detectable. 






















































10 




85 73 60 45 23 85 73 60 45 23 85 73 60 45 23 



FIG. 5: Amplitudes detectable at 3cr using temperature-only (left), polarization-only [center) and combined temperature and 
polarization [right) data from WMAP7 [orange), Planck [red) and a cosmic-variance-limited experiment [dark red). Detectable 
values of 7lo are indicated with solid lines, IZq with long-dashed lines and TIq with short-dashed lines. 


the factor of appearing in the linear amplitude, means the uncertainties in both parameters grow to the point 

where order-of-magnitude estimates are no longer possible, even with a full-sky, infinite-resolution, cosmic-variance- 
limited experiment. 


C. Comparing Datasets and Experiments 

e ; coll 

As discussed above, the detectability boundaries correspond to constant-amplitude curves in the —(1 — 

cos Axsep) plane. These amplitudes are plotted as a function of Xc in Fig. [^for all experiments considered. Fig. 
contains three panels, displaying the amplitudes detectable using temperature-only, polarization-only and combined 
temperature and polarization information. 

Concentrating first on the overall trends, we note that when only temperature information is used the detectable 
amplitudes decrease with Xc, but when only polarization information is used the detectable amplitudes in general 
increase with Xc- This is because the primordial curvature templates are normalized to the same central value, but 
the same is not true for the evolved CMB templates: there is slight variation in the central amplitude of the CMB 
templates [Lf^,Qf^) with Xc- In particular, the maximum values of Lj^ and increase with Xc, whereas the 
maximum values of and decrease with Xc- This behavior can be explained from the shape of the real space 
transfer functions used in the template computation (Eq. on large scales. In Fig. we display the functional 
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xlO-2 xlO"® 



r [Mpc] 

FIG. 6: At large scales, the temperature and polarization transfer functions behave qualitatively differently. While for polariza¬ 
tion, they remain positive over the full width of the last scattering surface (black solid line, scale shown on the right y-axis), for 
temperature, a change of sign occurs, leading to partial cancellation of the curvature perturbation contribution to the template 
signal (blue solid line, scale shown on the left y-axis). 


form of the integrand r‘^ai(r) for the multipole moment f = 20, corresponding to angular scales of ^ 10°. While the 
transfer function is larger than zero over the full width of the last scattering surface in polarization, in temperature, 
it changes sign from being positive at the location of the peak of the photon visibility function to being negative 
at smaller radii. For strictly non-negative curvature perturbations, contributions to the CMB temperature template 
will therefore partially cancel out, generating larger template amplitudes for more quickly decaying perturbations 
(i.e., larger Xc). For polarization, on the other hand, no such cancellation occurs and the CMB template amplitude 
increases for decreasing 

Turning to TZq , we again note that the quadratic signature is only detectable for a small range of the Xc values 
considered. The detectable TZq is smallest (and the range of Xc over which TZq is detectable is largest) when combining 
temperature and polarization. The CMB temperature provides stronger constraints than the polarization. Considering 
IZq instead, the sharp polarization feature present at Oc in the linear template means that IZq, in contrast to TZq , 
is well constrained by polarization alone; indeed it is better constrained by polarization in the range Xc < 3 Gpc. 
Temperature is, however, much better for the total amplitude IZq (red), as the quadratic template is so poorly 
constrained by polarization. 

For the overall amplitude, IZq, the detectability limits are mostly determined by temperature data. Note that 
the linear amplitude is always better constrained than the overall amplitude when employing only polarization, but 
the opposite is true when using temperature-only data. This is because the amplitude of the linear polarization 
template (specifically, the sharp feature at Oc) is around twice that of the quadratic polarization template, whereas 
the temperature templates are very similar in amplitude (Qj^ is typically ^ 10% stronger than Lj^). The IZo^f 
degeneracy is therefore exposed at larger 1 — cos Axgep when using polarization than when using only temperature. 
When temperature and polarization data are combined, the detectable IZq and IZq curves are very similar. 

Comparing experiments, the WMAP, Planck and cosmic-variance-limited curves all converge when using 
temperature-only information: both WMAP and Planck are essentially cosmic-variance-limited in temperature at 
the multipoles relevant to this analysis. Indeed, Planck is also nearly cosmic-variance-limited in E polarization, 
yielding constraints on TIq almost indistinguishable from the cosmic-variance-limited experiment, and constraints on 
IZq only 7% weaker at the largest angular scales. Though PlancPs polarization noise becomes more important for 
smaller collision signatures, its combined temperature and polarization data yield detectable amplitudes a maximum 
of only 8% larger than the cosmic variance limit on the smallest scales considered here. As the constraints on IZq 
are temperature-dominated, Planck offers less than a 10% improvement over WMAP; however, its polarization data 
allow up to ^ 30% improvement in constraints on the linear amplitude, IZq, for collision signatures covering half the 
sky. Beyond increasing the fraction of foreground-cleaned sky available for analysis, there is little improvement to be 
gained from post-Planck data. 
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V. COMPARING WITH PREVIOUS WORK 

The forecasts presented above are broadly consistent with previous work on constraining bubble collisions using data 
from the WMAP satellite. Two groups have presented constraints on the linear portion of the template. The most 
recent analysis of Feeney et. al. |45] implemented a hierarchical Bayesian formalism which tests the bubble collision 
model on the full sky, constraining the expected total number of collisions to be fewer than 4 at 95% confidence. To 
facilitate comparison with the single-template forecast presented above, we consider their candidate-detection step 
which implements an optimal filtering of the CMB for templates of different scales. The threshold for 50% of bubble 
collisions to be detected with a signal-to-noise ratio of greater than 3 is AT/T{uj) > 6.5 x 10“^, and is relatively 
independent of the angular scale of the collision Oc- Neglecting the evolution of the comoving curvature perturbation, 
in the Sachs-Wolfe approximation, where AT/T 7^/5, this implies IZq > 3.3 x 10“^. We take this to be the rough 
constraint on individual templates implied by the results of Feeney et. al. Osborne et. al. m implemented an 
optimal estimator for single collision templates, directly constraining IZq to |7^o I 10“"^(1 — cos^c)(sin^c)~^^^ 

at 2(7.^ Both of these constraints are compatible with the WMAP forecast presented in Fig. however the presence 
of sky cuts and other complications lead to somewhat weaker constraints from WMAP data than the full-sky forecast 
presented above. 

A competing method for detecting bubble collisions is to use the technique of kinetic Sunyaev-Zekdovich tomog¬ 
raphy [56]. This method circumvents cosmic variance by looking at the effect of a collision on the velocity field of 
free electrons, measured using the cross-correlation between large scale structure and the kinetic Sunyaev-Zekdovich 
effect in the CMB. A combination of Large Synoptic Survey Telescope (LSST) [57] and Planck temperature data is 
forecasted to yield constraints of IZq < 2 x 10“^ and 71q < 7.5 x 10“^ at 3cr for collisions on the largest angular 
scales. These constraints can be improved by an order of magnitude by a future high-resolution CMB experiment, 
with an ultimate cosmic-variance limit of < 10“^. 


VI. DISCUSSION AND CONCLUSIONS 

We have forecasted the ability of cosmic-variance-limited CMB datasets, employing information from both temper¬ 
ature and polarization, to constrain theories of eternal inflation using bubble collisions. We have studied both the 
raw detectability for bubble collisions, as well as the ability of CMB datasets to constrain fundamental parameters of 
the theory underlying eternal inflation. Comparing the detectability achievable with a cosmic-variance-limited CMB 
experiment and data from the Planck satellite to existing constraints arising from WMAP temperature data, we find 
that temperature data yields virtually no improvement on current limits. Adding polarization data improves limits on 
detectability by < 30%, and there are only marginal improvements between Planck and the cosmic-variance-limited 
experiment. We can therefore view data from Planck as having the final word on bubble collisions in the CMB. 

The bubble collision template for the comoving curvature perturbation in Eq. once evolved into the CMB temper¬ 
ature and polarization anisotropies, directly links the scalar field theory underlying eternal inflation to cosmological 
observables. In this paper, we have restricted ourselves to collisions between identical bubbles, in which case the 
observables are determined entirely by the width of the barrier in the scalar-field potential (fixed for a given theory) 
as well as the initial separation between bubbles and the observer position (stochastic variables). Assuming the de¬ 
tection of negative spatial curvature and primordial gravitational waves, which enter into the amplitudes for different 
components of the template, we have presented limits on detecting the barrier width, initial separation, and observer 
position. In the optimistic scenario where negative spatial curvature and primordial gravitational waves are detected, 
we expect these to be close to the bounds achievable by an analysis of Planck data. 

Looking beyond the CMB, it is possible to use measurements of the kinetic Sunyaev-Zekdovich effect to mine 
cosmological data further for signatures of bubble collisions. From the theoretical side, more-varied phenomenology 
could arise in theories of eternal inflation with multiple scalar fields, or for collisions between different bubbles. 
For example, recent work [58ll6Q] has shown a variety of models will break the assumed symmetry of the collision 
spacetime, possibly leading to a collision signature with more structure. We leave an exhaustive study of constraints 
on other models to future work. 


This is obtained by relating the parameter a in Ref. [47] to by = a(a^is - Xc). 
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